THEORY OF e'/e 



STEFANO BERTOLINI 
INFN and SISSA, Via Beirut 4, 1-34013 Trieste, Italy 
E-mail: bertolin@he.sissa.it 

I shortly review the present status of the theoretical estimates of e'/e . I consider a 
few aspects of the theoretical calculations which may be relevant in understanding 
the present experimental results. In particular, I discuss the role of final state 
interactions and in general of non-factorizable contributions for the explanation 
of the A J = 1/2 selection rule in kaon decays and e'/e . Lacking reliable "first 
principle" calculations, phenomenological approaches may help in understanding 
correlations among theoretical effects and the experimental data. The same dy- 
namics which underlies the CP conserving selection rule drives e'/e in the range 
of the recent experimental measurements. 

The results announced this year by the KTeV0 and NA48I collaborations 
have marked a great experimental achievement, establishing 35 years after the 
discovery of CP violation in the neutral kaon systemu the existence of a much 
smaller violation acting directly in the decays. 

While the Standard Model (SM) of strong and electroweak interactions 
provides an economical and elegant understanding of indirect (e) and di- 
rect (e') CP violation in term of a single phase, the detailed calculation of 
the size of these effects implies mastering strong interactions at a scale where 
perturbative methods break down. In addition, CP violation in K — ► irn 
decays is the result of a destructive interference between two sets of contri- 
butions, which may inflate up to an order of magnitude the uncertainties on 
the hadronic matrix elements of the effective four-quark operators. All that 
makes predicting e'/e a complex and subtle theoretical challengeu. 

The status of the theoretical predictions and experimental data available 
before the KTeV announcement early this year is summarized in Fig. ^. 

The gray horizontal bands show the one-sigma average of the old (early 
90's) NA31 (CERN) and E731 (Fermilab) results. The vertical lines show the 
ranges of the published theoretical predictions (before February 1999), iden- 
tified with the cities where most of the group members reside. The range of 
the naive Vacuum Saturation Approximation (VSA) is shown for comparison. 

The experimental and theoretical scenarios have changed substantially 
after the KTeV announcement and the subsequent NA48 result. Fig. shows 
the present experimental world average for e'/e compared with the revised or 
new theoretical calculations that appeared during this year. 

Notwithstanding the complexity of the problem, all theoretical calcula- 
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Figure 1. The 1-cr results of the NA31 and E731 Collaborations (early 90's) are shown by 
the gray horizontal bands. The old Miinchen, Roma and Trieste theoretical predictions for 
e' /e are depicted by the vertical bars with their central values. For comparison, the VSA 
estimate is shown using two renormalization schemes. 



tions show a remarkable overall agreement, most of them pointing to a non- 
vanishing positive effect in the SM (which is by itself far from trivial) . 

On the other hand, if we focus our attention on the central values, many 
of the predictions prefer the 10 -4 regime, whereas only a few of them stand 
above 10~ 3 . Is this just a "noise" in the theoretical calculations? 

The answer is no. Without entering the details of the various estimates, 
it is possible to explain most of the abovementioned difference in terms of a 
single effect: the different size of the hadronic matrix element of the gluonic 
penguin Qq obtained in the various approaches. 

While some of the calculations, as the early Miinchen and Rome predic- 
tions, assume for (irir\Qs\K) values in the neighboroud of the leading 1/N 
result (naive factorization), other approaches, among which the Trieste and 
Dortmund calculations, find a substantial enhancement of this matrix element 
with respect to the simple factorization result. The bulk of such an effect is 
actually a global enhancement of the 1 = components of the K — > tttt am- 
plitudes, which affects both current-current and penguin operators, and it can 
be clearly understood in terms of chiral dynamics (final-state interactions). 

As a matter of fact, one should in general expect an enhancement of 
e'/e with respect to the naive VSA due to final-state interactions (FSI). As 
Prof. Suzuki reminded usO, in two body decays the isospin 1 = final states 
feel an attractive interaction, of a sign opposite to that of the 1 = 2 compo- 
nents. This feature is at the root of the enhancement of the 1 = amplitude 
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Figure 2. The recent theoretical calculations of e'/e are compared with the combined 1-a 
average of the NA31, E731, KTeV and NA48 results (e'/e = 21.2 ± 4.6 X 10~ 4 ), depicted 
by the gray horizontal band (the error is inflated according to the Particle Data Group 
procedure when averaging over data with substantially different central values). 



over the 1 — 2 one and of the corresponding enhancement of e'/e beyond 
factorization. 

The question is how to make of a qualitative statement a quantitative one, 
since a reliable calculation of e'/e must involve at the same time a convincing 
explanation of the A7 =1/2 selection rule, which is still missing. n 

The AI = 1/2 selection rule in K — > tttt decays is known since 45 yeardS 
and it states the experimental evidence that kaons are 400 times more likely 
to decay in the 7 = two-pion state than in the 7 = 2 component. This 
rule is not justified by any general symmetry consideration and, although it is 
common understanding that its explanation must be rooted in the dynamics 
of strong interactions, there is no up to date derivation of this effect from first 
principle QCD. 

The ratio of I = 2 over 7 = amplitudes appears directly in the definition 

of e'/e : 

e' _ 1 \ ({™) I=2 \H W \K L ) ((™)i=2\Hw\K s ) \ 

e V2 \((™)i=o\Hw\K L ) ((Tnr)i=o\Hw\Ks) J ' U 

As a consequence, a consistent calculation of e'/e must also address the de- 
termination of the CP conserving amplitudes. 

The way we approach the calculation of the hadronic K — ► ttk transi- 
tions in gauge theories is provided by the Operator Product Expansion which 
allows us to write the relevant amplitudes in terms of the hadronic matrix 
elements of effective AS = 1 four quark operators (at a scale /i) and of the 
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corresponding Wilson coefficients, which encode the information about the 
dynamical degrees of freedom heavier than the chosen renormalization scale: 

G 

72 



Has=i = -j= v ud v* s + T Vi(v) Qi(jJ>) 



(2) 



The entries Vij of the 3x3 Cabibbo-Kobayashi-Maskawa matrix describe the 
flavour mixing in the SM and r = —VtdV t */V u dV* s . For fi < m c (q — u, d, s), 
the relevant quark operators are: 



Qi = (s a u p ) v _ A {u f} d a ) Y _ A 1 
Q 2 = {su) v _ A (ud) v _ A J 

Qs,5 = (sd) v _ A Y, q (qq)v T A 

Qi,6 = (s Q d / 3) v _ A I] q (9/39a)vTA 



{7,9 = | (sgOv-a Eg e g (qq) 



V±A 



Current-Current 



Gluon "penguins" 



Electroweak "penguins'' 



(3) 



V±A 



Current-current operators are induced by tree-level W-exchange whereas the 
so-called penguin (and "box" ) diagrams are generated via an electroweak loop. 
Only the latter "feel" all three quark families via the virtual quark exchange 
and are therefore sensitive to the weak CP phase. Current-current operators 
control instead the CP conserving transitions. This fact suggests already 
that the connection between e'/e and the AI =1/2 rule is by no means a 
straightforward one. 

Using the effective AS — 1 quark Hamiltonian we can write e'/e as 



where 



e 

e 

n 
n 2 



j<t>_ 



GfLlI 



2|e|ReA 



Im A f 



n 



1 



n 2 



cos <5q 



cos 8' A 



(4) 



(5) 



and (Qi) = (inr\Qi\K). The-rescattering phases 60,2 can be extracted from 
clastic 7r — 7r scattering dataQ and are such that cos^o ~ 0.8 and cos(52 — 1. 
Given that the phase of e, 9 S , is approximately 7r/4, as well as 5q — 62, <j> = 
■| + 82 — 5q — 9 e turns out to be consistent with zero. 
Two key ingredients appear in eq. |]: 

1. The isospin breaking ir° — rj — 77' mixing, parametrized by f2^+^', which 
is estimated to give a positive correction to the A2 amplitude of about 
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15-35%. However, it has been recently emphasized!! that NLO chiral cor- 
rections may make VL v+ri > negative and as large as —0.7, thus potentially 
providing a strong enhancement mechanism for e'/e . 

2. The combination of CKM elements ImAt = Im(V t * s Vtd), which affects 
directly the size of e'/e and the range of the uncertainty. .The determi- 
nation of Im At goes through B-physics constraints and e The latter 
observable depends on the theoretical determination of Bk , the K° — K a 
hadronic parameter, which should be consistently determined within ev- 
ery approach. The theoretical uncertainty on Bk affects subtantially the 
final uncertainty on Im At . A better determination of the unitarity trian- 
gle is expected from the B-factories and the hadronic collidersEJ. From 
K-physics, Kl — > ir°vD gives the cleanest "theoretical" determination of 
Im At , albeit representing a great experimental challenged. 

A satisfactory approach to the calculation of e'/e should comply with the 
following requirements: 

A: A consistent definition of renormalized operators leading to the correct 
scheme and scale matching with the short-distance perturbative analysis. 

B: A self-contained calculation of all relevant hadronic matrix elements 
(including Bk)- 

C: A simultaneous explanation of the AI =1/2 selection rule and e'/e . 

None of the available calculations satisfies all previous requirements. I 
summarize very briefly the various attempts to calculate e'/e which have 
appeared so far and have lead to the estimates shown in Figs. |l| and |2|. A 
simple naive approach to the problem is the VSA, which is based on two 
drastic assumptions: the factorization of the four quark operators and the 
vacuum saturation of the completeness of the intermediate states. As an 
example: 

(ir+7r-\Q 6 \K°) = 2 (tt" \uj 5 d\0) (ir+ \su\K°) - 2 (Tr + ir-\dd\0)(0\s l5 d\K°) 

+ 2 [(0|ss|0) - (0\dd\Q)] (ir+ir~\sj 5 d\K°) (6) 

The VSA does not allow for a consistent matching of the scale and scheme 
dependence of the Wilson coefficients (the HV and NDR results are shown in 
Fig. |l|) and it carries a potentially large systematic uncertainty!! It is best 
used for LO estimates. 

Generalized Factorization represents an attempt to address the issue of 
renormalization scheme dependence in the framework of factorization by defin- 
ing effective Wilson coefficients which absorb the perturbative QCD running 
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of the quark operators. The new coefficients, scale and scheme independent at 
a given order in the strong coupling expansion, are matched to the factorized 
matrix elements at the scale /ip which is arbitrarily chosen in the pertur- 
bative regime. A residual scale dependence remains in the penguin matrix 
elements via the quark mass. Fitting the AI — 1/2 rule and e'/e requires 
non-factorizablej-contributions both in the current-current and the penguin 
matrix elementsEJ. 

In the Munchen approach (phenomenological 1/N) some of the matrix 
elements are obtained by fitting the AI = 1/2 rule at /j, = m c = 1.3 GeV. 
On the other hand, the relevant gluonic and electroweak penguin (Qq) and 
{Qb) 2 remain undetermined and are taken around their leading 1/N values 
(which implies a scheme dependent result). In Fig. || the HV (left) and NDR 
(right) results are shownEl The dark range represents the result of gaus- 
sian treatment of the input parameters compared to flat scannning (complete 
range) . 

In the recent years the Dortmund group has revived and improved the ap- 
proach of Bardeen, Buras and GerarcO based on the 1/N expansion. Chiral 
loops are regularized via a cutoff and the amplitudes are arranged in a p 2n /N 
expansion. A particular attention has been given to the matching procedure 
between the scale dependence of the chiral loops and that arising from the 
short-distance analysistH The renormalization scheme dependence remains 
and it is included in the final uncertainty. The AI =1/2 rule is reproduced, 
but the presence of the quadratic cutoff induces a matching scale instability 
(which is very large for Bk)- The NLO corrections to (Qe) induce a substan- 
tial enhancement of the matrix element (right range in Fig. ||) compared to 
the leading order result (left). The dark range is drawn for central values of 
m s , ImA( and A QC d- p 

In— .the Nambu, Jona-Lasinio (NJL) modelling of QCDE3 the Dubna 
groupliH has calculated e'/e including chiral loops up to 0(p 6 ) and the effects 
of scalar, vector and axial-vector resonances. Chiral loops are regularized via 
the heat-kernel method, which leaves unsolved the problem of the renormal- 
ization scheme dependence. A phenomenological fit of the AI = 1/2 rule 
implies deviations up to a factor two on the calculated {Q§). The reduced 
(dark) range in Fig. |^ corresponds to taking the central values of the NLO 
chiral couplings and varying the short-distance parameters. The A J = 1/2 
rule and Bk have been, studied in the NJL framework and 1/N expansion 
by Bijnens and Pradealj showing an impressive scale stability when including 
vector and axial-vector resonances. 

In the approach of the Trieste group, based on the Chiral Quark Model 
(xQM)E3, all hadronic matrix elements are computed up to 0(p 4 ) in the chi- 
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ral expansion in terms of the three model parameters: the constituent quark 
mass, the quark condensate and the gluon condensate-irhese parameters are 
phcnomenologically fixed by fitting the AI = 1/2 rulec3. This step is crucial 
in order to make the model predictive, since there is no a-priori argument 
for the consistency of the matching procedure. As a matter of fact, all com- 
puted observables turn out to be very weakly dependent on the scale (and the 
renormalization scheme) in a few hundred MeV range around the matching 
scale, which is taken at 0.8 GeV as a compromise between the ranges of valid- 
ity of perturbation theory and chiral expansion. The 1 = matrix elements 
are strongly .enhanced by non-factorizable contributions and drive e'/e in the 
10 -3 regimeEil. 

Lattice regularization of QCD is the consistent approach to the prob- 
lem. On the other hand, there are presently important numerical and the- 
oretical limitations, like the quenching approximation and chiral symmetry, 
which may substantially affect the calculation of the weak matrix elements. 
In addition, chiral perturbation theory is needed in order to obtain K — ► irn 
amplitudes from the computed K — > tt transitions. As summarized by Soni 
at this conferencecj lattice cannot provide us at present with reliable calcula- 
tions of the 1 = penguin operators relevant to e'/e , as well as of the 1 = 
components of the hadronic matrix elements of the current-current operators 
(penguin contractions), which are relevant to the AI = 1/2 rule. This is due 
to large renormalization uncertainties, partly related to the breaking of chiral 
symmetry on the lattice. In this respect, very promising is the Domain Wall 
Fermion approachc 3 i-which allows us to decouple the chiral symmetry from 
the continuum limit£3. In the recent Romapte-evaluation of e'/e (Qe) is taken 
at the VSA value with a 100% uncertainty^.! The result is therefore scheme 
dependent (the HV and NDR results are shown in Fig. ||). The dark (light) 
ranges correspond to Gaussian (flat) scan of the input parameters. 

Other attempts to reproduce the measured e'/e using the linear cr-model, 
which include the effect of a scalar resonance with m a ~ 900 MeV, obtain 
the needed enhancement of (Qq)cQ. However, it is not possible to reproduce 
simultaneously the experimental values of e'/e and of the CP conserving K — > 
7T7T amplitudes. 

Without entering into the details of the various calculations I wish to 
illustrate with a simple exercise the crucial role of final-state interactions (in 
general of non-factorizable contributions) for e'/e and the AI = 1/2 selection 
rule. In order to do that I focus on two semi-phenomenological approaches. 

A commonly used way of comparing the estimates of hadronic matrix 
elements in different approaches is via the so-called B factors which represent 
the ratio of the model matrix elements to the corresponding VSA values. 
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Figure 3. Predicting e'/e: a (Penguin) Comparative Anatomy of the Miinchen (dark gray) 
and Trieste (light gray) results (in units of 10 — 3 ). 



However, care must be taken in the comparison of different models due to 
the scale dependence of the -B's and the values used by different groups for 
the parameters that enter the VSA expressions. An alternative pictorial and 
synthetic way of analyzing different outcomes for e'/e is shown in Fig. ||, where 
a "comparative anatomy" of the Trieste and Miinchen estimates is presented. 

From the inspection of the various contributions it is apparent that the 
different outcome on the central value of e'/e is almost entirely due to the 
difference in the Qe component. 

In the Munich approacfJld the AI —1/2 rule is used in order to determine 
phenomenologically the matrix elements of Qi,2 and, via operatorial relations, 
some of the matrix elements of the left-handed penguins. The approach does 
not allow for a phenomenological determination of the matrix elements of 
the penguin operators which are most relevant for e'/e , namely the gluonic 
penguin Q 6 and the electroweak penguin Q 8 . These matrix elements are taken 
around their leading 1/N values (factorization). 

In the semi-phenomenological approach of the Trieste group the size of 
the effects on the 7 = 0,2 amplitudes is controlled by the phenomenological 
embedding of the AI = 1/2 selection rule which determines the ranges of 
the model paremeters: the constituent quark mass, the quark and the gluon 
condensates. In terms of these parameters all matrix elements are computed. 

Fig. |i| shows an anatomy of the xQM contributions which lead to the 
experimental value of the AI = 1/2 selection rule for central values of the 
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Figure 4. Anatomy of the AI = 1/2 rule in the xQM. See the text for explanations. The 
cross-hairs indicate the experimental point. 



input parameters. 

Point (1) represents the result obtained by neglecting QCD and taking the 
factorized matrix element for the tree- level operator Q 2 , which is the leading 
electroweak contribution. The ratio A j A 2 is thus found equal to V2: by far 
off the experimental point (8). Step (2) includes the effects of perturbative 
QCD renormalization on the operators 2t£l. Step (3) shows the,effect of 
including the gluonic penguin operatorscll. Electroweak penguinaHS are nu- 
merically negligeable in the CP conserving amplitudes and are responsible for 
the very small shift in the A 2 direction. Therefore, perturbative QCD and 
factorization lead us from (1) to (4). 

Non-factorizable gluon-condensate corrections, a crucial model dependent 
effect entering at the leading order in the chiral expansion, produce a substan- 
tial reduction of the A 2 amplitude (5), as it was first observed by Pich and 
de RafaelnJ. Moving the analysis to 0(p 4 ), the chiral loop corrections, com- 
puted on the LO chiral lagrangian via dimensional regularization and minimal 
subtraction, lead us from (5) to (6), while the finite parts of the NLO coun- 
terterms calculated in the xQM approach lead to the point Cty. Finally, step 
(8) represents the inclusion of ir-rj-r]' isospin breaking effectsEI 

This model dependent anatomy shows the relevance of non-factorizable 
contributions and higher-order chiral corrections. The suggestion that chiral 
dynamics may be relevant to the understanding of the A T —1/2 selection rule 
goes back to the work of Bardeen, Buras and GeraroO in the 1/N framework 
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Figure 5. Anatomy of e'/e (in units of 10 — 3 ) in the xQM approach. In black the LO 
results (which includes the non-factorizable gluonic corrections), in half-tone the effect of 
the inclusion of chiral-loop corrections and in light gray the complete 0(p 4 ) estimate. 



with a cutoff regularization. A pattern similar to that shown in Fig. [f| for the 
chiral loop corrections to A$ and Ai was previously obtained in a NLO chiral 
lagrangian_analysis, using dimensional regularization, by Missimer, Kambor 
and WylerEJ. The Trieste group has extended their calculation jXLinclude the 
NLO contributions to the electroweak penguin matrix element s! 3 ^rH 

Fig. H shows the contributions to e'/e of the various penguin operators, 
providing us with a finer anatomy of the NLO chiral corrections. It is clear 
that chiral-loop dynamics plays a subleading role in the electroweak penguin 
sector (Qs—w) while enhancing by 60% the gluonic penguin (7 = 0) matrix 
elements. The NLO enhancement of the Qq matrix element is what drives 
e'/e in the xQM to the 10~ 3 ballpark. 

As a consequence, the xQM analysis shows that the same dynamics that 
is relevant to the reproduction of the CP conserving Aq amplitude (Fig. ||) is 
at work in the CP violating sector, albeit with a reduced strenght. 

In order to ascertain whether these model features representj-ceal QCD 
effects we must wait for future improvements in lattice calculations^. On the 
other hand, indications for such a dynamics follow from the analitic properties 
of the K — > 7T7T amplitudes. We should expect in general an enhancement of 
e'/e , with respect to the naive VSA estimate, due to final-state interactions. 
In two body decays, the 1 = final states feel an attractive interaction, of a 
sign opposite to that of the 1 = 2 components. This feature is at the root of 
the enhancement of the 1 = amplitude over the 1 = 2 one. Recent dispersive 
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analyses^ of K — > tttt amplitudes show how a (partial) resummation of final 
state interactions increases substantially the size of the 1 = components, 
while slightly depleting the 1 = 2 components. 

It is important to notice however that the size of the effect so derived 
is generally not enough to fully account for the AI = 1/2 rule. Other non- 
factorizable contributions are needed to further enhance the CP conserving 
1 = amplitude, and to reduce the large 1 = 2 amplitude obtained from 
perturbative QCD and factorization. In the xQM approach, for instance, the 
fit of the AI = 1/2 rule is due to the interplay of FSI (at the NLO) and 
non-factorizable soft gluonic corrections (at LO in the chiral expansion) . 

The idea of a connection between the AI = 1/2 selection rule and 
e'/e goes back a long waya, although at the GeV scale, where we can trust 
perturbative QCD, penguins are far from providing the dominant contribution 
to the CP conserving amplitudes. 

I conclude by summarizing the relevant remarks: 

Those semi-phenomenological approaches which reproduce the AI =1/2 
selection rule in K — > 7T7t decays, generally agree in the pattern and size of 
the 1 = 2 hadronic matrix elements with the existing lattice calculations. 

1 = amplitudes: the AI =1/2 rule forces upon us large deviations from 
the naive VSA: B— factors of O(10) for (Q\^)q- Lattice calculations presently 
suffer from large sistematic uncertainties. 

In the Trieste and Dortmund calculations, which reproduce the CP con- 
serving K — > 7T7T amplitudes, non-factorizable effects (mainly final-state inter- 
actions) enhance the hadronic matrix elements of the gluonic penguins, and 
give Bq/B^P « 2. Similar indications stem from other l/iV@ and dispersiveE! 
approaches. Promising work is in progress on the lattice. 

Theoretical error: a substantial amount of the present uncertainty in the 
e'/e prediction arises from the uncertainty in the CKM elements lm(V t *Vtd) 
which is presently controlled by the AS = 2 parameter Bk- A better determi- 
nation of the unitarity triangle from B-physics is expected from the B-factories 
and hadronic colliders. From K-physics, =r ir°vD gives the cleanest "the- 
oretical" determination of Im X t . 



New Physics: in spite of recent clever proposals (mainly SUSYE3) it is 
premature to invoke physics beyond the SM in order to fit e'/e . A number 
of ungauged systematic uncertainties affect presently all theoretical estimates, 
and, most of all, every attempt to reproduce e'/e must also address the puzzle 
of the AI =1/2 rule, which is hardly affected by short-distance physics. 
An "anomalously" large e'/e is likely to be the CP violating counterpart of 




A /A 2 » 22. 
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